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L'invention concerne un proceed de localisation d'un ou de 
5 plusieurs emetteurs. 

II s'applique notamment a la localisation d'emetteurs au sol a partir 
d'un engin mobile sans avoir de connaissance a priori sur les signaux emls. 

Le domaine technique est notamment celui de la localisation 

10 passive d'emetteurs. 

La figure 1 illustre une localisation aeroportee. L'emetteur est a la 
position (x 0l yo,z 0 ). Le porteur a I'instant t k est a la position (x*,y*,z k ) et voit 
l'emetteur sous ['incidence (6<fc,x 0 ,y 0 ,Zo) I A(t k ,xo,y 0 ,Zo)). Les angles 
6(t,xo,yo,zo) et A(t,x 0 ,yo,Zo) evoluent au cours du temps et dependent de la 

15 position de l'emetteur ainsi que de la trajectoire du porteur. 

Les angles ^f,x 0> y 0 ,z 0 ) et 4(f,x 0 ,y 0> z 0 ) sont reperes par rappbrt a 
un reseau de N antennes pouvant etre fixe sous le porteur comme le montre 
la figure 2. 5 

II existe actuellement plusieurs techniques permettanf; de 

20 determiner la position (x m , y m , z m ) d'un emetteur. Ces techniques' de 
localisation different notamment par les parametres qui sont estimes en 
instantanee au niveau du reseau des capteurs. Elles peuvent etre classees 
de la maniere suivante: 
Utilisation de la goniometrie, 

25 Ces techniques sont connues et utilisees dans I'art anterieur. 

Dans la plupart des cas, elles sont basees sur une goniometrie 1D en 
azimut. Les azimuts & n =6{ffoX m ,y, nf z m ) associes au m iSme emetteur sont 
mesures pour differents instants f fc En utilisant la position (x k ,y k ,z k ) du porteur 
a I'instant correspondant k, une position (x^.y^,^) de l'emetteur m est 

30 estimee par une intersection sol. La position (xioVk.zK) du porteur est donnee 
par un GPS, son orientation est obtenue par un compas dans le cas d'un 
porteur terrestre et par une centrale de navigation dans le cas d'un aeronef. 
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A partir de toutes les positions (x mk ,y mk ,z mk ), la methode effectue une 
extraction de donnees permettant de determiner les M positions dominantes 
(x m ,y m ,z m ) des emetteurs incidents. La localisation est obtenue par 
triangulation ou par intersection sol (goniometrie 2D). L'inconvenient des 
5 techniques de triangulation est qu'elles necessitent un defilement important. 
D'autre part, les techniques de goniometrie doivent utiliser un reseau de 
capteurs non ambigue pour fournir les incidences. Ceci a pour inconvenient 
de necessiter une table de calibration et de limiter la taille du reseau de 
capteurs et par consequent de fournir des incidences limitees en precision. 
10 Utilisation de la difference de phase entre 2 capteurs eloignes, 

La difference de phase A<p(t k ,Xo,yo,Zo) entre capteurs depend de la 
positions des 2 capteurs ainsi que de I'incidence (<%,x 0 ,y 0 ,z 0 ), d(fc,x 0 ,y 0 ,zo)) 
de I'emetteur. Cette phase qui depend du temps est directement liee a la 
position (xo,y 0 ,z 0 ) de I'emetteur. En consequence, en etudiant la fonction du 
15 temps A<p(t,xo,yo,zo) il est possible d'en deduire la position (x 0 ,y 0 ,z 0 ) de 
I'emetteur. Dans cette famille d'application les 2 capteurs sont eloignes pour 
augmenter la precision de la mesure de la phase. Ceci a pour inconvenient 
de faire varier la difference de phase ^£xo,yo,z 0 ) en fonction du temps sur 
plus de 2n et la technique necessite alors une etape permettant de derouler 
20 la phase sur plus de 2rc. D'autre part dans cette technique la phase est 
mesuree en effectuant directement une intercorrelation entre 2 capteurs ce 
qui ne permet pas de traiter le cas multi-emetteurs. 
Utilisation de la mesure de la frequence porteuse de I'emetteur, 

Ces techniques exploitent le fait que la frequence porteuse 
25 estimee est la somme de la frequence porteuse de I'emetteur et du decalage 
doppler du a la vitesse de deplacement du porteur. Le decalage doppler a 
I'avantage de dependre de la position (x 0 ,y 0 ,z 0 ) de I'emetteur et d'etre aussi 
une fonction du temps Af(f,x 0 ,yo,z 0 ). En consequence en etudiant la fonction 
du temps Af(f,x 0 ,y 0l z 0 ) il est possible d'en deduire la position (x 0 ,y 0 ,z 0 ) de 
30 I'emetteur. La mesure de ce decalage doppler presente toutefois comme 
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inconvenient de necessiter des emetteurs ayant des formes d'ondes 
particulieres. Cette mesure de frequence peut se faire par des techniques 
cycliques supposant que Ie signal emis est non circulaire. 
Utilisation des temps de propagation, 
5 Ces techniques exploitent les differences de temps de propagation 

entre aeriens (TDOA ou Time difference of arrival) qui sont directement liees 
aux distances respectives de I'emetteur aux differents aeriens et done a la 
position (Xo,y 0 ,Zo) de l'6metteur. En utilisant au moins trois aeriens 
suffisamment espaces, il est possible de deduire ia position (x 0l yo,Zo) de 

10 I'emetteur par localisation hyperbolique. ^inconvenient de ces techniques est 
qu'elles ne peuvent etre mises en oeuvre en contexte mono porteur en raison 
des espacements considerables requis entre aeriens. D'autre part dans r ces 
techniques, la difference de temps est mesuree en effectuant directement 
une intercorrelation entre 2 capteurs, ce qui ne permet pas de traiter Iq< cas 

15 multi-emetteurs. 

.?* 

Le precede selon Tinvention repose notamment sur *une 
nouvelle approche d'estimation directe des positions (x m , y m , z m ) de chapun 
des emetteurs a partir d'une analyse parametrique du signal multi-voies a 
20 divers instants tk sur une duree Dt. L'analyse parametrique a notamment 
pour fonction supplementaire de separer les differents emetteurs a chaque 
instant tk- On associe ensuite les parametres d'un meme emetteur issus des 
differents instants tk pour finalement localiser chacun des emetteurs. 

25 ^invention concerne un precede de localisation d'une ou de 

plusieurs sources, la ou lesdites sources 6tant en mouvement par rapport a 
un reseau de capteurs, le precede comportant une etape de separation des 
sources afin d'identifier les vecteurs directeurs associes a la reponse des 
capteurs a une source d'incidence donnee. II est caracterise en ce qu'il 

30 comporte au moins les etapes suivantes : 
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• associer les vecteurs directeurs a^.-.a,^ du m'eme emetteur obtenus 
respectivement aux instants t 1 ...t K> 

• localiser le m i6me emetteur a partir des vecteurs a., ...a„ associes. 

1 m Km 

Le precede selon invention presente notamment les avantages 

suivants : 

• il permet de localiser en plus de la position en (x, y, z) d*un emetteur son 
vecteur vitesse, 

• il s'applique lorsque I'on est en presence de un ou plusieurs emetteurs 
incidents, 

• sa mise en oeuvre ne necessite pas de connaissances particulieres sur le 
signal emis, 

• il permet d'utiliser un reseau de capteurs ambigue c'est-a-dire plusieurs 
incidences sont associees a la meme reponse du reseau qui ont 
I'avantage d'etre grand et ainsi d'etre plus robuste aux phenomenes de 
couplage entre a<§riens ou plus generalement aux erreurs de modelisation 
du reseau d'aeriens, 

• il peut etre mis en oeuvre sur des reseaux calibres en (9,A). 

• il peut etre mis en oeuvre sur des reseaux a antennes a diversites 
d'amplitude comme les antennes colocalisees : reseau avec des dip6les 
de meme centre de phase et ayant des orientations differentes. 

D'autres caracteristiques et avantages de I'objet de la presente 
invention apparaTtront mieux a la lecture de la description qui suit donnee a 
titre illustratif et nullement limitatif a la lecture des figures annexees qui 
representent : 

• la figure 1 le schema de principe de la localisation d'un emetteur de 
position au sol au moyen d'un aeronef, 

o la figure 2 la relation entre un reseau d'antennes et I'incidence d'un 
6metteur, 



• la figure 3 un schema general expliquant le fonctionnement du precede 
selon I'invention, 

• les figures 4, 5 et 6 des exemples de mise en oeuvre du procede selon 
I'invention. 

Afin de mieux faire comprendre I'objet de la presente invention, la 
description qui suit est donnee a titre illustratif et nullement limitatif pour 
iocaliser plusieurs emetteurs disposes au sol au moyen d'un reseau de 
capteurs equipant un aeronef en mouvement. Un tel systeme est par 
exemple decrit a la figure 1. L'aeronef est equipe d'un processeur adapte a 
mettre en ceuvre les etapes du procede selon I'invention. 

Le procede peut aussi §tre mis en oeuvre dans le cadre de 
vehicule en mouvement au sol. . 

La figure 3 represente, dans un diagramme temps-amplitude .du 
signal, le signal x(t) compose d'une combinaison des signaux des emetteurs 
a differents instants t,, t 2 , ...t K . Sur cette figure les differentes etapes mises 
en ceuvre a savoir la separation des emetteurs SE et I'estimation 
parametrique EP, I'association des parametres de chaque emetteurji la 
localisation d'un emetteur sont resumees. 

En presence de M emetteurs, le procede dispose, a I'instant t en 
sortie des N capteurs du reseau, du vecteur x(f) representatif du melange 
des signaux des M emetteurs. Autour de I'instant t k , le vecteur x(f+f*) de 
dimension Nxt, representant le melange des signaux des M emetteurs, 
s'exprime de la maniere suh/ante : 

M ... 
*(t + tk) = 2 a (©ton. A/cm) s m (t+t k )+ b(t+t k )=A k s(t+t k ) + b(t+t k ) pour \t\<At/2 ( } 

ou b(f) est le vecteur bruit suppose gaussien, a(8,A) est la reponse du reseau 
de capteurs a une source d'incidence (9,A), A^fafe/r,, A fc1 )... 3(9^, Aw)], 
s(0=[si(0--SM(f)3 T . ^i< m =0it k ,x m ,y m ,z m ) et A km =A(t k ,x m ,y m ,z m ). Dans ce modele 
la matrice A^ de melange depend de I'instant t k d'observation. 
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Le vecteur directeur de I'incidence correspondant au /n i6me emetteur a 
I'instant t k 

a km = a(6 km , A /(n ,)=a(f fo x m ,y mi z m ) du m i6me emetteur (2) 

est une fonction connue de t k et de la position de I'emetteur (x m ,y m ,z m ). 

5 Le procede selon ('invention comporte au moins les etapes 

suivantes : 

1 . estimer un ou plusieurs parametres associes a la position de la source, 
par exemple les vecteurs directeurs, les incidences, la position, etc. et 
separer les M emetteurs pour les differents instants t^, ce qui consiste a 

10 identifier les vecteurs directeurs d'incidence a km pour (1<m<M). Cette 
premiere etape est par exemple effectuee par des techniques de 
separation de sources connues de PHomme du metier, 

2. associer les parametres estimes pour le m'^me emetteur, par exemple en 
associant les differents vecteurs directeurs d'incidences, a-j m a km 

15 obtenus respectivement aux instants t-f,...^, 

3. localiser le m i6me emetteur a partir des vecteurs associes. 
Etape dissociation 

En presence de M emetteurs et apres separation de sources, le 
procede possede a I'instant t k les M signatures a^ pour (1<m<M). A I'instant 
20 t k +, la separation de source donne les M vecteurs b, pour (1</<M)- L'objectif 
de ce suivi est de determiner pour le m i6me emetteur, I'indice i(m) qui 
minimise I'ecart entre a^ et b i(m) . Dans ce cas on en deduira que a k+1> m = 
bi (m) . Pour effectuer cette association on definit par exemple la distance entre 
deux vecteurs u et v par : 



w ^ < l U " V \ ^ 

(a H «J[v H v ) 

25 OD u H est le transpose conjugue du vecteur u. 
Dans ces conditions I'indice i(m) verifie : 
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d(a*™ b m ) = min [d^, b/)] (4) 



done 



H. 



d(a* m , = min[1- £ i 7\ J 



Dans cette association on considere une fonction p m associee au m i6me 
emetteur: 

Pm(tk) = d(a km , aom) (5) 

5 Au fil de I'association on obtient pour chaque emetteur m et pour -\<m<M, la 
fonction p m (t). Cette fonction a notamment pour objectif d'eliminer les instants 
t k dont la valeur p^fc) paraTt trop eloignee d'une interpolation de la fonction 
p m (f), e'est-a-dire que Ton elimine les instants aberrants qui peuvent §tre 
associes a d'autres emetteurs. On definit une zone de tolerance +/-A agtour 
10 de la courbe definie par la fonction $ m (t k ). Cette zone de tolerance dependra 
de la precision d'estimation des vecteurs directeurs. a*™. En particulier en 
presence de M=1 source la zone sera de I'ordre de A=3/^/5At (ou At est le 
temps elementaire d'estimation parametrique illustre fig.3 et B est la bande 
instantanee du signal x(t) ). 

15 Les etapes de cette association pour K instants t k sont par 

exemple les suivantes : 

Etape ASE - 1 : Initialisation du processus a k=2. Le nombre M d'emetteurs 
initial est par exemple determine par un test de detection du nombre de 
sources a I'instant t 0 connu de I'Homme du metier, 
20 Etape ASE - 2 : Pour 1<m<M determination des indices i(m) en appliquant 
I'equation (4) et en utilisant le vecteur a Km avec Um<M et les vecteurs b, 
identifies a I'instant f* +1 pour (1<i<M), 

Etape ASE - 3 : Pour 1 <m<M effectuer I'operation a k+1 m = b i(m)i 
Etape ASE - 4 : Incrementation k<-k+1 et si k<K retour a I'etape ASE-1 , 
25 Etape ASE - 5 : A partir de la famille d'instants <D={ f,<. . .< t K }, eliminer les / 



8 



instants f,e<t> tel que les coefficients n'appartiennent pas a la zone 
delimitee par la courbe d'interpolation des Mt k ) et fa zone de tolerance A. 
On elimlnera aussi les instants tk ou IMtk)- Pm(t k -i)|<A. Apres ce tri la 
nouvelle famille d'instants est 0={ U<...< t, } et on pose K=l. 
5 A la fin de ces etapes, le precede a determine les vecteurs a 1/n 

a/c m associes au m ,6me emetteur. 

Localisation d'un emetteur 

Le precede determine la position du ni 6 ™ emetteur a partir des 
composantes des vecteurs a 1m jusqu'a a Km . Ces vecteurs a Am ont la 
10 particularite de dependre de I'instant t k et surtout de la position (x m ,y m ,z m ) de 
I'emetteur. En particulier pour un reseau compose de N=2 capteurs espaces 
d'une distance de tfdans I'axe du porteur le vecteur verifie a*™: 

a =f ( d 1 \] (6) 

La valeur 1 de la premiere composante correspond au capteur de 

reference. D'apres la figure 1, I'incidence (m^,y m ,z m ), A{t^ m ,y mi 7i m ) ) peut 
15 etre directement calculee a partir de la position (x k ,y k ,z k ) du porteur a I'instant 

t k et la position (Xm.y^z^ de I'emetteur. 

Etape de transformation du vecteur 

Selon une variante de realisation, le precede comporte une etape 

de transformation du vecteur a^ en un vecteur b^ dont les composantes 
20 sont formees a partir des composantes du vecteur a^. En particulier, le 

precede construit par exemple le vecteur b^ de dimension (A/-1)x1 en 

choisissant un capteur de reference en n=i : 
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bkm = 



»j-(*-i)/«*.(0 

(' + !)/«*. (0 



(7) 



b(4,x m ,y m ,z m ) ou a^C/) est la / 6me 



composante de a*™ 

Les composantes de b* m correspondent dans ce cas aux rapports des 

composantes du vecteur et du vecteur a^i). 

Ainsi dans I'exemple de i'equation (6) en fixant /=1 on obtient : 



exp\j2xj.cos(0(t k , Xm ,y m , Zm >co S (A( ?i ,x m ,y m , z m >J 



(8) 



= a Am (2)/a ton (1) 

Sachant que les vecteurs directeurs a km sont estimes avec une 
certaine erreur ekm tel que a* m = a(f to x m ,y m ,z m )+ efen , on peut en deduire.qu'il 
en est de meme pour le vecteur transforme b^, de (7) . ,> : 
Etape de maximisation d'un critere de correla tion 

Sachant que le vecteur a km est une fonction de la position 
(Xm,ym,Zm) de I'emetteur il en est de meme pour le vecteur b km . Le procede 
comporte une etape de maximisation d'un critere de correlation vectorielle 
normalise Ux.y.z) dans I'espace (x.y.z) de position d'un emetteur ou 

Ux,y,zK r * w (9) 



Avec 



b K = 



'Im 



_ Km _ 



= VKCXm.ym.Zn,) + W K , V/c (X.y.Z) = 



10 



Ira 



et wk = 



w 



Km 



Le vecteur bruit w K a pour matrice de covariance R=E[w K w K H ]. En 
faisant I'hypothese que la matrice R est connue, le critere peut etre envisage 
avec une technique de blanchiment. Dans ces conditions on obtient le critere 
Lk'(x,v,z) suivant : 



Avec R=E[w k Wk h ] 

5 » II faut remarquer que les criteres des equations (9) et (10) sont egaux 
lorsque R=a 2 l, c'est a dire lorsque les erreurs sont considerees de niveau 
egal sur tous les capteurs et ind^pendantes entre capteurs. Le critere 
LK(x,y,z) de I'equation (10) est done valable pour un bruit de modele w K de 
. statistiques blanches. 

10 • Les criteres L^x.y.z) et L K '(x,y,z) sont compris entre 0 et 1 et verifient 
L K (x,y,z)= L K '(x,y,z)=1 pour la position (Xm.ym.Zn,) du /77 idme emetteur. Cette 
normalisation permet de fixer un seuil de bonne localisation t\. Ainsi tous les 
maximums (Xm.ym.Zm) de L^x.y.z) qui verifient L^Xm.ym.Zn,^ t) sont 
consid§res comme des bonnes localisations. Le seuil peut etre fixe en 

15 fonction d'une connaissance approchee des statistiques de w K . 

« Les criteres L/<(x,y,z) et Lx'foy.z) ont I'avantage de pouvoir mettre en 
oeuvre une technique de localisation en presence d'un reseau de capteurs 
calibres dans I'espace (9,A). Sachant qu'a I'instant t k on connaTt la relation 
analytique liant I'incidence (e(f/ f> x,y,z), A(^,x,y,z)) de I'emetteur a sa position 

2p (x,y,z), on peut alors deduire a partir de I'incidence (e^.x.y.z), A(^,x,y,z)) le 
vecteur a(f /( ,x n , 1 y mi z m )= a(Q(t k ,x,y,z), A(t k ,x,y,z)) en realisant une interpolation 
de la table de calibration (relative aux antennes calibrees). Remarquons 



L K '(x I y,z)= 



(10) 



^ H R"'b,>, ix t y,z)*R-*Y K (x,y, z) ) 
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cependant que ce procede est insensible a un biais en phase (en raison du 
critere de correlation vectorielle). 

• Ces criteres permettent aussi de tenir compte de la phase et de I'amplitude 
des composantes de a(0,A). La rhethode peut done §tre envisagee avec des 
5 reseaux a antennes colocalisees a diversite de diagramme. 

II faut remarquer que dans un contexte aeroporte la connaissance de 
Paltitude h de I'avion permet de reduire le calcul du critere dans I'espace de 
recherche (x,y) en posant z=/?. Dans I'exemple des equations (6) et (8) le 
vecteur v/<(x,y,z) s'ecrit de la maniere suivante : 

(11) 



v K (x,y,z) = 



exp|y2?r-|-cos(0(f , , x, y, z)) cos(A(*, , x, y, z)) j 
exp|y27z- j cos(B(t K , x, y, z> cos( A(^ , x, y, z> j 



10 Dans ce procede il est possible d'envisager d'initialiser 

ralgorithme a K=K 0 puis ensuite de calculer de facon recursive le critere 
L/<(x,y,z). Dans ces conditions L/c(x,y,z) se calcule recursivement de la fagon 
suivante : » 



L K .,<x,y, Z )= K-<*** 



(12) 



ou 

aK + i(x,y,z)= a^x.y.z) + b K+im H b(f K+ i,x,y,z) 
YK + i(x,y,z)=YK(x,y,z) + b(f K+ i,x,y,z) H b(fe + i,x f y f z) 

Pk+1 =Pk + bK+1 m H b K +i m 

Les coefficients a K+ i(x,y,z)= a,<(x,y,z) y K+1 (x,y f z)=y K (x ) y,z) ) p K+ i=P* 
15 sont des spectres intermediaires permettant de calculer l_K+i(x,y,z). 

Lorsque les vecteurs b(f/c + i,x,y,z) e t b^ sont de normes 
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constantes egale a p\a relation de recurrence de I'equation (12) devient : 
ou 

aK + i(x,y,z)= a K (x,y,z) + b K+ * m H b(f K+1 ,x,y,z) 

Le procede est decrit jusqu'ici en supposant que les emetteurs ont 
des positions fixes. II peut facilement s'etendre au cas de cibles mobiles de 
vecteur vitesse (v^Vy^Vzn,) pour lesquelles on dispose d'un modele 
devolution. Dans ces conditions I'incidence du m i6me emetteur se 
parametrise de la maniere suivante : 

Qkm=8(t k , X^-Vxm t k ,y m -Vym t k , Z^zm t k ) et A krri =A( t k , X m -V xm t k ,y m -V ym t k , Z„rV zm t k ) 

ou (x m ,y m ,z m ) est la position de I'emetteur a I'instant t 0 et (Vxm.Vym.Vzm) les 
composantes de la vitesse de I'emetteur a I'instant t 0 . Dans ces conditions le 
vecteur de I'equation (7) est parametre par (w^z*,) et (v^.v^.v^) de 
la maniere suivante : 

b/ox, = b(f fo x m ,y m ,z m ,v xm ,Vym,v zm ,)+ w frm 

De maniere naturelle les criteres de localisation L K et L K ' des 
equations (9) et (10) ne sont plus parametres seulement par (x,y,z) mais 
aussi par (v Xl v y ,v z ). Le procede consiste done a maximiser le critere L,<(x,y,z, 
v x ,Vy,v z ) en fonctions des 6 parametres (x.y.z, v x ,v y ,v z ). 

Le procede peut s'appliquer a un tres grand nombre de mesures. 
Dans ce cas, le procede comporte une etape de reduction de la complexity 
numerique de calcul (qui est fonction du nombre de mesures) en diminuant 
K. Le procede prevoit d'effectuer sur les mesures elementaires les 
traitements suivants : 

- decimation des instants t k , en eliminant les instants voisins pour lesquels 
revolution de la courbe p m (f fc ) n'est pas significative, 

- filtrage (lissage des mesures qui sont les vecteurs directeurs) et sous 



13 



echantillonnage, 

- les mesures sont ensuite fusionnees sur une duree definie (extraction par 
association de vecteur directeur pour produire une mesure de synthese). 

Recapitulation des etapes du procede 

5 Le procede de localisation de plusieurs emetteurs utilisant K instants t k 

peut se resumer par les etapes suivantes : 

Etape n°1 : Identification des vecteurs a km pour (1 <m<M) aux K instants f k en 
appliquant par exemple une technique de separation de sources et 
d'identification de sources comme decrit dans les references [2] [3]. 
10 Etape n°2 : Association des vecteurs a im jusqu'a a Km obtenus aux instants 
respectifs ft... t K associe au m l6me emetteur pour 1<m<M en appliquant les 
etapes ASE-1 jusqu'a ASE-5 decrit ci-dessus. 
Etape n°3 : Initialisation du processus a /7?=1 

Etape n°4 : Transformation des K vecteurs a km en des vecteurs b K m comme 

. x . • • f 

15 le suggere Pequation (7). 

Etape n°5 : Calcul et maximisation du critere L.K(x,y,z) de Tequation (9) pour 

obtenir la position (x™,y m ,z m ) du rri*™* emetteur. 

Etape n°6 : Incrementation m<-/n+1 et si m<M retour a Tetape n°3 

Afin d'affiner I'estimation de la position (x m ,y m ,y m ) des emetteurs les etapes 
20 du procede peuvent etre realisees de fagon iterative de la maniere suivante : 

Etape n°7 : Identification des vecteurs b, pour (1</</W) a Tinstants f K +i en 

appliquant par exemple une technique de separation et d'identification de 

sources comme decrit dans les references [2] [3]. 

Etape n°8 : Pour 1<m<M determination des indices i(m) en appliquant 
25 I'equation (4) et en utilisant le vecteur a*™ et les vecteurs b/ pour (1 <i<M). 
Etape n°9 : Pour ^m<M est effectue Poperation a K +i m= b i(m) 
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Etape n°10 : Pour 1<m<M calcul du critere U +1 (x,y,z) de facon iterative en 
utilisant les equations (12) et (13) et minimisation de L K+ i(x,y,z) pour obtenir 
la position (Xm^Zm) du m i6me emetteur. 

Etape n°11 : Si Ton decide de continuer pour etre plus precis et moins 
ambigue le procede retourne a I'etape n°7. 

Exemple de mise en oeuvre du procede 

Les simulations ont ete realisees avec un reseau de N=2 capteurs 
alignes dans I'axe du porteur avec d/X=3. Comme dA=3 une methode 
effectuant une goniometrie aux instants t k serait completement ambigue et ne 
permettrait pas de faire par la suite des triangulations pour effectuer la 
localisation de I'emetteur. Sur les figures 5, 6 et 7 correspondant aux criteres 
de localisation pour K=3, 7 et 16 est trace le pseudo-spectre L K (x,y) a 
maximiser permettant de determiner la position de I'emetteur dans I'espace 
(x,y). Sachant que si I'emetteur se situe en (x 0 ,y 0 ) alors L K (xo,yo)=1, on en 
ded.uit que les courbes iso-niveaux L K (x,y)=0.99 caracterisent la largeur du 
lobe principale. Remarquant que la precision de localisation depend de la 
largeur de ce lobe, on en deduit d'apres les figures 5, 6 et 7 que plus K est 
important et meilleur sera la precision de localisation. 
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REVENDICATIONS 



1 - Procede de localisation d'un ou de plusieurs sources, la ou lesdites 
sources etant en mouvement par rapport a un reseau de capteurs, le 

5 procede comportant une etape de separation des sources afin d'identifier les 
vecteurs directeurs associes a la reponse des capteurs a une source 
d'incidence donnee, caracterise en ce qu'il comporte au moins les etapes 
suivantes : 

• associer les vecteurs directeurs a 1m -..a KlTl obtenus pour le m^me 
10 emetteur et respectivement pour les instants ^...t^,, 

• localiser le mi© 1118 emetteur a partir des vecteurs a 1m ...a Km associes. 

2 - Procede selon la revendication 1 caracterise en ce que Tetape 
dissociation comporte au moins les etapes suivantes : 

15 Etape ASE - 1 : initialiser le processus a k=2, 

Etape ASE - 2 : pour A<m<M determiner les indices i(m) en utilisant la 
relation ,d(a/cm, b i(m) ) = nunfdCa^ £>/)], le vecteur a ktm et les vecteuts b, 

identifies a Tinstant t k ^ pour (1</<M), etablir une fonction p m (t k ) = d(a km , a&m) 
Etape ASE - 3 : pour 1 <m<M effectuer i'operation a^ m= b,y m ;, 
20 Etape ASE - 4 : incrementer /c«-/c+1 et si k<K retourner a Tetape ASE-1 , 
Etape ASE - 5: a partir de la famille d'instants 0={ fi<...< t K }, ainsi 
obtenue, extraire les instants f/ qui n'appartiennent pas a une zone definie 
par la courbe p m (tk) et une zone de tolerance. 

25 3 - Procede selon la revendication 1 caracterise en ce que Tetape de 
localisation comporte au moins les etapes suivantes : 

maximiser un critere de correlation vectorielle normalise L k (x, y, z) dans 
Tespace (x, y, z) de position d'un emetteur avec 
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REVENDICATIONS 



jl - Procede de localisation d'une ou de plusieurs sources, la ou lesdites 
sources etant en mouvement par rapport a un reseau de capteurs, le 
procede comportant une etape de separation des sources afin d'identifier les 
vecteurs directeurs associes a la reponse des capteurs a une source 
d'incidence donnee, caracterise en ce qu'il comporte au moins les etapes 
suivantes : 

• associer les vecteurs directeurs a 1m ...a Km obtenus pour le m ieme 
emetteur et respectivement pour les instants t...t 

1 K* 

• localiserlemieme Emetteur a partir des vecteurs a, ...a, associes. 

ifn Km 

2 - Procede selon la revendication 1 caracterise en ce que I'etape 
d'association comporte au moins les etapes suivantes : 

Etape ASE - 1 : initialiser le processus a k=2, 

Etape ASE - 2 : pour 1<m<M determiner les indices i(m) en utilisant la 
relation d(a^ b l(m) ) = nuB[d( aton , *>,)], le vecteur a kifn et les vecteurs b, 

identifies a I'instant t k+ , pour (1 <i<M), etablir une fonction p m (tk) = d(a km , a om ) 
Etape ASE - 3 : pour 1<m<M effectuer I'operation a k+1m = b i(m)t 
Etape ASE - 4 : incrementer k<r-k+1 et si k<K retour a I'etape ASE-1 , 
Etape ASE - 5 : a partir de la famille d'instants <D={ ^<...< t K }, ainsi 
qbtenue, extraire les instants t, qui n'appartiennent pas a une zone definie 
par la courbe p m (tk) et une zone de tolerance. 

3 - Procede selon la revendication 1 caracterise en ce que I'etape de 
localisation comporte au moins les etapes suivantes : 

maximiser un critere de correlation vectorielle normalise Lk(x, y, z) dans 
I'espace (x, y, z) de position d'un emetteur avec 
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Avec 



b K = 



L b *»r 



= v/cCx^y^Zm) + w K , v K (x,y,z) = 



b(t l9 x 9 y,z)' 
b(t K9 x 9 y y z) 



et wk - 



w 



Km 



oil w/c est le vecteur bruit pour toutes les positions (x, y, z) d'un emetteur. 

4 - Procede selon la revendication 3 caracterise en ce que le vecteur b K 
comporte un vecteur representatif du bruit dont !es composantes sont 
fonctions des composantes des vecteurs ai m . . . ai<m. 



10 5 - Proced6 selon la revendication 3 caracterise en ce qu'il comporte une 
etape ou I'on determine la matrice de covariance R=E[wk w k h ] du vecteur 
bruit et en ce que Ton maximise le critere 



U (x,y,z)= -r— - — ; — ^ — J * 



6 - Procede selon la revendication 5 caracterise en ce que revaluation du 
15 critere L/^x^z) et/ou du critere L.K(x,y,z) est recursive. 

7 - Proced6 selon Tune des revendications 1 a 6 caracterise en ce qu'il 
comporte une etape de comparaison des maximums avec une valeur seuil. 

20 8 - Procede selon Tune des revendications 1 a 7 caracterise en ce que la 
valeur de K est initialement fixee a K 0 . 
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K M v *(*,J'..*) H v*(;c,;,,z) 



5 



Avec 



- VK(x m ,y m ,z m ) + w K , v K (x,y,z) = 



et Wk = 



w 



in 



w, 



6u w* est le vecteur bruit pour toutes les positions (x, y, z) d'un emetteur. 

4 - Procede selon la revendication 3 caracterise en ce que le vecteur b K 
comporte un vecteur represented du bruit dont les composantes sent 
fonctions des composantes des vecteurs a 1m ... aKm. 

5 - Procede selon la revendication 3 caracterise en ce qu'il comporte une 
etape ou Ton determine la matrice de covariance R=E[w K w K H ] du vecteur 
bruit et en ce que Ton maximise le critere 

6 - Procede selon la revendication 5 caracterise en ce que revaluation du 
critere Ufx.y.z) et/ou du critere L K '(x,y,z) est recursive. 



7 - Procede selon I'une des revendications 1 a 6 caracterise en ce qu'il 
comporte une etape de comparaison des maximums avec une valeur seuil. 
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9 - Precede selon Tune des revendications 1 a 8 caracterise en ce que les 
emetteurs a localiser sont mobiles et en ce que le vecteur considere est 
parametre par la position de Pemetteur a localiser et ie vecteur vitesse. 
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8 - Precede selon I'une des revendications 1 a 7 caracterise en ce que les 
emetteurs a localiser sont mobiles et en ce que le vecteur considere est 
parametre par la position de I'emetteur a localiser et le vecteur Vitesse. 
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